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1 . Perform statistical analyses and correlation studies, and develop classification 
schemes, e g., histograms, two-color diagrams, correlation plots, classification of spectral 
features 

2. Construct phenomenological models for each class of objects to parameterize the 
observed rartafon characteristics to terms of their physical source properties. 

3. From the phenomenological models, develop self-consistent physical models which can 
explain coherently present and future observations. 



flNPUTl i OUTPUT 


• luminosity and spectral energy distribution 
of heat source 

• source geometry & dimensions 

• properties of grain species (e g., opacity, 
albedo, scattering phase function) 

• density distribution of dust grains 

• other local heating mechanisms (e.g., 
viscosity, collisions with gas particles) 


* temperature distributions of each grain 
species 

* characteristics of internal radiation field 

* intensity variation across source surface at 
different wavelengths 

* characteristics of grain spectral features 

* emergent energy spectrum of heat source 
and dust shell 


OBSERVATIONAL CONSTRAINTS & DIAGNOSTIC PROBES 


• spectral energy distribution (Buk spectrum) o source luminosity 

• surface brightness & apparent emission size as a function of wavelength » dust density 
distribution & geometrical source size 

• strength and shape ot emission/absorption features o grain composition & source opacity 

• color temperatures m different parts of spectrum o dust temperature distribution 

• flux ratios or colot -color diagrams o source evolution 


Figure 1 - Schematic diagram to show that a crucial step in the theoretical 
studies of infrared sources is the construction of phenomenological models. 


To model the observed characicrislics of an infrared source, typically one solves 
the equation of radiation transport in a dusly medium, subject to the conslraini of 
radiative equilibrium For spherically symmetric geometry, the problem has been 
solved by many authors (e g , Leung 1975, Apruzesc 1976; Jones and Merrill 1976; 
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2. SOME RECENT ADVANCES 
2. 1 Grain Properties 

In aslrophysical environments dust grains have irregular shapes most likely formed by 
fractal growth processes, e.g., particle-by-particle aggregation or cluster-by-cluster 
(Witten and Cates 1986). On the other hand, spherical dust grains are 
often assumed in models of infrared sources so that the dust opacity can be calculated 
from the Mie theory. A computational technique now exists for calculating the dust 
opacity for grains of irregular shapes (Draine 1988; Bazell and Dwek 1990). In this 
method, called "discrete dipole approximation", a fractal grain is approximated by a 
collection of dipoles. Using realistic dust opacities for fractal grains, recently Fogel 
and Leung (1992) have constructed radiation transport models for infrared sources. 
Compared to models with spherical grains of the same composition and volume, 
models with fractal grains show a shift in the peak flux toward longer wavelengths, 
implying that fractal dust grains are cooler than spherical grains . These differences 
can be attributed to a lower ratio (p) of volume to geometric cross section (averaged 
over orientation) for the less compact fractal grains. The ratio p varies with the fractal 
dimension. Grains with a higher p attain a higher temperature. Furthermore, grains 
having the same fractal dimension show almost no difference in their absorption 
cross sections and energy spectrum, implying that the overall grain shape plays only 
a minor role in the thermal properties . 

In modeling infrared sources, the opacity of dust grains is often assumed to be 
temperature-independent. However, the optical constants of many grain types actually 
change with grain temperature. In particular, for water ice, which has a feature at 3.1 
^m, the absorption coefficients change with temperature. The ice feature is seen in the 
spectra of molecular clouds and circumstellar dust shells of both young and evolved 
stars (Whittet 1992). As the temperature increases, water ice changes from an 
amorphous state to a crystalline state. This irreversible phase change, which occurs at 
around 100 K, narrows the ice feature and increases the peak strength, thus producing 
different profiles for the 3.1 pm feature as the temperature increases. Since the grain 
temperature changes with position in infrared sources, it is crucial that temperature- 
dependent effects in grain opacity be incorporated in detailed models to interpret 
observations of ice features. This will allow us to probe the evolution of ice mantle on 
grains. Recently laboratory data on temperature-dependent optical constants for water 
icc (Hudgins cl al 1992) have become available, making it possible to include these 
effects in radiation transport models. Work is underway to study this problem. 

2.2 Physical Processes 

Among the uncxpecied results from IRAS was the discovery of excess mid-infrared 
emission detected in many infrared sources, e g., diffuse clouds, dark globules, visual 
reflection nebulae, and high-latitude dust clouds or infrared cirrus. It is now believed 
that the emission at short wavelengths (< 30 pm) comes from transient heating of very 
small grains and large polycyclic aromatic hydrocarbons or PAHs (for a review see 
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Figure 3 - Schematic diagram comparing the theoretical and observational 
considerations in modeling infrared sources with 1-D and 2-D geometries. 


2.4 Modeling Approach 

Unlike ground-based observations, space missions generally yield large quantities of 
data Interpreting such data requires a different strategy which utilises extensive 
modeling. As an example, the IRAS survey has yielded a large and uniform sample of 
high-quality data from observations of asymptotic giant branch stars. In particular, 
the IRAS data indicates that many carbon stars, especially those with optically thin 
dust shells, have large fluxes at 60 and 100 pm. It has been suggested that cither an 
extended single dust shell or a two-shell model, with a remnant shell from an earlier 
mass-loss episode (Willems 1987), can explain the excess fluxes. To test these 
hypotheses, models of dust around carbon stars have been constructed: models with 
cither a single C-rich dust shell or double shells (Egan and Leung 1991). Figure 4 
shows the color-color diagrams for selected carbon stars in the IRAS-LRS catalog. 
The thick solid line is the blackbody line. The gray-shaded areas indicate the limits on 
the colors for carbon stars of various ages imposed by the models. Clearly single-shell 
models cannot explain the observed color distribution, while two-shell models with 
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Pugel and Leger 1989). Temperature fluctuations in small grains occur whenever the 
energy inpul from photons or energetic particles is considerably larger than the 
average energy content of the grain. This non-equilibrium grain heating can change 
significantly the energy distribution of the radiation field in infrared sources. By 
formulating the radiation transport problem involving transient heating in a fashion 
similar to a non-LTE line transfer problem involving many transitions, Lis & Leung 
(1991a, b) recently constructed models which treat self-consistently the thermal 
coupling between the transient heating of small grains and the equilibrium heating of 
conventional large grains. They used the models to interpret the IRAS observations of 
the Barnard 5 cloud and a diffuse cloud in Charaaeleon. In both cases, longward of 
100 pm, the emission is dominated by large grains. Between 30-100 pm, the emission 
is produced mainly by very small grains. Shortward of 30 pm, both PAHs and small 
grains are responsible for the emission. Typically very small grains and PAHs 
account for 10-20% of the total opacity in the visible, and 5-20% of the total dust 
mass of the cloud. Furthermore, to produce the observed infrared limb brightening 
the spatial distribution of small grains and PAHs must be more extended than that of 
large grains. This has important significance in understanding the origin of small 
grains and PAHs. Thus detailed radiation transport models now exist which 
incorporate both the transient healing of very small grains and the equilibrium heating 
of conventional large grains- 

Another important progress deals with the study of grain formation in stellar 
oullfows, using a master-equation approach (Egan and Leung, this volume). This 
method allows a self-consistent treatment of grain nucleation and growth. It is found 
that, compared to the classical nucleation theory, the new method predicts fewer but 
larger grains . This has important observational implications. 

2.3 Source Geometry 

Although there is growing observational and theoretical evidence for a large number of 
disk-shaped or toroidal objects of aslrophysical interest (e.g., circumstellar disks, 
proioplanetary disks, proiostcllar accretion disks, bipolar molecular flows, and disk 
galaxies), most models curremly in use invoke the assumption of spherical geometry 
A few attempts have been made to model infrared sources with nonspherical geometry 
Lefevre el al. (1983) performed Monte Carlo simulations of ellipsoidal dust shells 
around cool siars, while Ghosh & Tandon (1985) calculated dust temperature 
distributions in cylindrical clouds with embedded stars. The latter work has been 
extended by Dent (1988) to the case of circumstellar disks around young stars. Most 
recently, Spagna et al (1991) considered radiation transport in disk-shaped interstellar 
clouds heated externally by the ambient interstellar radiation field. While theoretical 
tools are available for modeling infrared sources with different geometries, the typical 
thousand fold increase in computing requirement for 2-D geometries makes it 
impractical to do any extensive modeling. In Figure 3 we compare the two geometries 
(1-D sphere vs. 2-D disk) from theoretical and observational considerations. Since 
establishing source geometry places severe constraints on the origin, dynamics, and 
properties of infrared sources, it is crucial to consider realistic geometries. 
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3 CONCLUDING REMARKS AND OUTLOOK 
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can model infrared sources on workstations in real time. When computer modeling 
becomes as easy as doing least-square fits, it will become a routine part of data 
analysis. With automation researchers can perform computer experiments to test 
various hypotheses and determine the physical parameters for infrared sources. With 
such a research tool, predictions on certain observational consequences can be made 
readily. This will stimulate further observations and theoretical studies, efforts which 
are essential to the future space-based missions such as ISO and SIRTF. 
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